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Abstract. We study the semiclassical Wigner-Kirkwood (WK) expansion of the 
^-j- ' partition function Z(t) for arbitrary even homogeneous potentials, starting from the 

C^l ■ Bloch equation. As is well known, the phase-space kernel of Z satisfies the so-called 

Uhlenbeck-Beth equation, which depends on the gradients of the potential. We perform 
a chain of transformations to obtain novel forms of this equation that invite analogies 
with various physical phenomena and formalisms, such as diffusion processes, the 
Fokker-Planck equation, and supersymmetric quantum mechanics. 
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1. Introduction 



The problem of the semiclassical expansion of the physical states and their energies 



in power of h rather than the expansions involving oscillating functions of Ti 1 has a 
long history and was studied and used in different fields of science (see [l] for review) 
be ginning from the Thomas- Fermi (TF) approximation jlJEllH]. Being the lowest order 
in h, the TF correction to the classical partition function (heat kernel) takes into account 
^ ■ the discreteness of the phase-space treating the Hamiltonian classically. 

We consider here Hamiltonians of the form H = T + V where T is the kinetic and 
V the potential energy, respectively. The non-commutativity of the T and V in the 
quantum Hamiltonian 

H = -y V 2 + V({a*}) (1) 

where V 2 denotes the Laplacian and % labels the coordinates, was taken into account 
by the Wigner-Kirkwood (WK) approach jUE]. (Note: We set the mass m = 1, and 
all quantities are given in units of energy E (implying the dimensions: [H] = 1, [t] = 
— 1, [h] = 3/4, [V] = l,[x] = 1/4).) We do not specify the form of the, generally non- 
central, potential V(x) for n degrees of freedom (x = {x\, x 2 , ■ ■ ■ , x n )). The case of even 
homogeneous polynomials satisfying V(Ax) = \ 2N V(x) includes the potentials of the 
well-known Yang-Mills classical and quantum mechanics [3 13 E] (see jHj for a review) . 

f Present address: 3106 Hornbuckle Place, Durham, NC 27707, USA. 
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2. The Partition Function in the WK expansion 

The partition function 

oo 

Z(t) = Tr [exp(-fiT)] = £ e~ Ent , (2) 

n=0 

which is the Laplace transform of the density of states C[p(E)], is calculated taking 
the Wigner transform of the quantum operator exp(-tH) and integrating it over the 
phase-space with measure 

dT = (27rh)- n d n xd n p, (3) 

where n is the dimension of the system. 

For the semiclassical expansion it is convenient to use the plane waves as a complete 

set: 

Z{t) = {2ixh)- n J dre~ ip x/h e~ t6 e ip x/h . (4) 

The WK method takes into account only part of the quantum features, but ignors the 
quantum fluctuations defined by the concrete form of the potential and dynamics (see 
[8j for the discussion of the role of the quantum fluctuations in Yang-Mills quantum 
mechanics). Following [TP] , we can account the non-commuting terms in WK approach 
setting 

e -tii e i P -x/h = e -tH e i P - x /h w ^ p . ^ = u ^ p . ^ 

where H is the classical Hamiltonian and the function W(x, p;t) is to be determined. 
The function u(x, p; t) satisfies the Bloch equation, the analogue of the time dependent 
Schrodinger equation: 

| + tf„ = (6) 

with the initial condition 

limu(x,p;0) = e ip - x/n , (7) 

corresponding to the initial condition W(x, p; 0) = 1. 

Boundary conditions on w(x, p; t) at \x\ = oo are in correspondence with the 
boundary conditions on W(x, p;i). From (jSJ) and © we obtain an exact equation 
for W(x,p;t), which we call the Uhlenbeck-Beth (UB) equation): 

dW h 2 



dt " 2 



h 

Next we expand W in powers of h 



v 2 - t(VT) + t 2 (vvy - 2t(vv) ■ v 

■ (V - tW)l W. (8) 



W = E (9) 

fc=0 
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and equate the terms with the same power of % on the both sides. The result is: 
dW k 1 r _ 



dt 



v 2 - t(VT) + r(w) 2 - 2t(w) ■ v 



w, 



k-2 



+ip ■ (V - tW)W k -!. 
with conditions Wo = 1, Wk = for k < 0. 



(10) 



3. New forms of the UB equation 

The expressions (jSJ) and (jlUJI can be written in more compact form, if we introduce the 
"covariant" derivative T> = V — t(VV): 

dW n 2 



dt 2 
or its recursive form 

dW k _ 1 r 
<9t 2 



r^ 2 2 * 






P + 7T P 


V 





{hv + i P ) 2 + P ^ 



ir. 



P 2 Vy fc _ 2 + 2*p ■ VW k ^ 



(11) 
(12) 



We note that the symbol p in (jllll2j) denotes a classical phase-space variable and not 
an operator. To our knowledge, the forms (JTTJ) and (JT2J) of eqs. (jHJ) and (fTUj) are novel 
and have not been given in the literature. 

In terms of the operator V, eq. (fTTj) resembles a Fokker-Planck equation (FPE) 
for W(x, p; t) with the diffusion constant h 2 and the constant drift vector — ihp. The 
relation to the FPE can be further elucidated by noting that the "vector potential" 
A = tW is a complete gradient and thus may be "gauged out" by the transformation 
W — ► exp(tV)W , yielding an alternative form of (fTTj) : 



dW h 2 



Ot 



V 2 + — p ■ V 



w - VW'. 



(13) 



If we interpret W'(x, p; t) as an one-time probability density and introduce the 
probability current (sometimes called the probability flux in the literature fH\ I12j): 



/i 2 / ^ 2 N 

j = -ih-pW - yvr = w f -inp - y v in w 



we may write (|T3jl in the form of the continuity equation 
dW 



dt 



+ v- J' 



-V(x)W'(x,p;t), 



(14) 



(15) 



where the potential term acts as a source term leading to a "probability loss" and 
violates the local conservation law associated with FPE. Integrating over volume and 
using Gauss' theorem in (|TH|) we obtain 

d_ 

dtJ js 

where S is the surface confining the volume. If the flux vanishes at infinity we have 
d 



J W'd n x + J s da.J' = -J VW'd 1 



x, 



(16) 



I W'd n x = - I VW'd n x, 



(17) 
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Note that in the right-hand side of (JT3J) and (fT5)l we may substitute the total energy H 
instead of the potential V by making the transformation W = W" exp(— p 2 t/2), since 
p commutes with gradients on the right-hand side of (12) and (14). 

The presence of the potential term VW (or UW) in ()13|) and ()15j) invalidates the 
standard result for the FPE (see, e. g. fUEj) that under rather loose conditions on 
the drift and diffusion kernels (which are satisfied in our case) all solutions of the FPE 
must coincide for sufficiently large times. If drift and diffusion functions do not depend 
on time, stationary solution may exist, and one solution is unique in the sense that all 
other solutions agree with it after sufficiently long time. Having the stationary solution 
of the FPE, one is able to find the general solution with given initial and boundary 
conditions (see, e. g. [HUE]). 

In our case it is clear that stationary solutions to the Bloch equation © exist only 
for the trivial case of vanishing V. This is evident from (jlj, © and (JBJ): Using the 
second expression for the probability flux in (JHJ), setting W st = and using U instead 
of V in (fTHjl . we obtain the equation 

(Vlnu st ) 2 = -V(x). (18) 
n 

Since u st = exp(zp • x/fo) is the solution of the Bloch equation © with du/dt = 0, 
we arrive at the stated conclusion that there are no stationary solutions to the Bloch 
equation for non-vanishing potential. Similarly, for the FPE (J13J) resulting from the 
substitution W = exp(tV)W , it is clear from (jSJ) that the only stationary solution is 
W = 1 for V(x) = 0. 

Kepping this in mind, let us divide the solution to (fT3j) into a stationary and a 
time-dependent part: 

W'(x, p; t) = <(x, p) + «/(x, p; t). (19) 

For V = there exists the stationary solution 



h 2 f_o 2i 



V 2 + -p ■ V 

n 



W' st = 0. (20) 



L FP (x,p)W/^ o 
Inserting (fl"9"j) into (fT^j) we obtain: 

^ = (L FP - V)w> - VW* (21) 

Following the standard approach to the solution of inhomogeneous first-order differential 
equations, we can obtain a formal solution of (|2~Tj) in the form 

u/(x,p;t) = - f dtV^-^-*V(x)^ t (x,p). (22) 

J — oo 

To this we need to add the general solution of the homogeneous equation 

^ = (L FP - V)w. (23) 

The complete solution ()19|) thus has three contributions: 

W'(x, p; t) = <(x, p) + «/(x, p; t) + £ C n w n (^ p; t), (24) 
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where the C n are derived from the initial condition W'(x, p; 0) = 1. The solutions 
w n of the homogeneous equation are obtained by separation of variables, u>„(x, p;t) = 
T n (t)X n (x, p), in the form 

T„(t) = e- A ™*, (25) 
(L - V0X„(x, p) = -A 2 X„(x, p), (26) 

where —A 2 are the eigenvalues of the operator L — V. Finally, the solution of (|2T)j) is 
the plane wave: 

l^ t (x,p) = ^e- 2 ^. (27) 
4. Chain of transformations of the UB equation 

By the transformation W = exp(tV)W (or W = exp(tU)W) we obtained the equation 
which is intermediate between the Bloch and UB equations. Going one step further 
by setting W = exp(— ip ■ sl/K)W we arrive, of course, at the Bloch equation (JHJ) as it 
is evident from the chain of transformations 

W -> e tu W" -> e 4l/ - ip - x/?i W, (28) 

with the identification W = u.) Each link of this chain opens up new analogies with 
the important physical processes and equations which govern them. We elaborate this 
issue below in more detail. Before we address this subject, however, we note that the 
expansion in h proposed by Kirkwood 5\ is also an expansion in powers of the gradient 
operator as emphasized in ref. ^U] : The gradient operator occurs k times for each Wk 
in eq. (|TUj) . In the general case, one needs to expand in powers of % or, equivalently, in 
powers of the gradient operator using (jTOjl and (fT2|) . as done in [7||H]. 

In some special cases the compact form of (fTTj) or (|13|1 and the Bloch equation 
itself may be the starting point of another approximation scheme, which can bring to 
bear the knowledge of approximation techniques developed for the Fokker-Planck and 
Schrodinger equations. With this in mind, we explore here several modified forms of the 
diffusion equation (|T3j). before we consider some aspects of the Bloch equation. Going 
back to our original function W(x, p; t) = exp(tU)W'(x, p; t) we may write (|T3|) in the 
form 

(>]V e tH (L FP W)e- tH , (29) 



dt 



where 



Lpp -y 



V 2 + — p ■ V 

h 



(30) 



For (JT5j) we then obtain 



dt 
with 



° ]V e tH (V ■ J)e~ tH , (31) 
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2 ' 2i 



W. (32) 



Equations and (|31|) are of form of a FPE, just like eq. (jSJ). 

Another, equivalent form of the same equation (fT3*j) is easily derived: 

^ = ie« - p*)e-« (33) 

with the quantum operator 

P = p - AV. (34) 

This is the Bloch-Schrodinger form of eq. (jSJ). 

As a check, from Z(t) = J dTWe~ tH we have, using (J3~3~j) : 

| = -/dr(^ + y(x))^e^. (35) 

Integrating over £ we have 

Z(t) = JdT {^P 2 + V{4) J°° dt'e-*' 11 W{x, P ; t'). (36) 

In the limit h = 0,V = p, (and thus = 1) eq. ()36|) reduces to the Thomas-Fermi 
term 

Z(t) = J dTe- tH . (37) 

Finally, we mention another form of the basic equation (jSJ), or its equivalent, 
eq. (fTTj) . If one introduces instead of W the new function W = exp(— tV) W, then 
(II lj) can be written in the form 

where the operator T> is given by 

V = e- tv / 2 Ve- ty l 2 . (39) 

Equation ()38|) corresponds to the FPE used in supersymmetric quantum mechanics in 
one dimension in the association with the Darboux transformation (see, e. g. 

5. More on the Bloch equation 

We see that WK semiclassical expansion leads to several equations with interesting 
connections to fundamental branches of physics. The basis of this expansion is the 
Bloch equation ©• Having obtained the solution of © we find the partition function 

Z(t) = J dr e - ip - x/n u(x,p;t). (40) 

The initial condition 

u(x, p; 0) = e ip - x/h (41) 

may appear somewhat unusual taking into account that the potentials V(x) considered 
here are generally associated with the bound states of the corresponding Schrodinger 



— V 2 + ihp-V 
2 P 



equation. Nevertheless, starting from we may represent w(x, p;t) as an expansion 
in terms of plane waves satisfying the initial condition (}4"T]) : 

u(x, p; t) = / dVa(p, p'; t) e * p '- x/?i - p ' 2 ' /2 . (42) 
Equation ()41j) translates into the following initial condition for the amplitude a(p, p'; t): 
a(p )P ';0) = 5(p-p'), (43) 

and eq. © gives, after multiplying the resulting equation by exp(— ip ■ x/ft) and 
integrating both sides over x: 



0_ 

Of 



a(p,p;t)e-^/ 2 ] = - / 7^a(p, p'; t)e~^ 

J d n xV{x)e i{p '- p) x/h . (44) 
Introducing the function 

A(p,p';t)=a(p,p';t)e-^/ 2 (45) 

we obtain 

-A(p, p; t) + P -A{p, p; t) = - J d n p'A(p, p', t)V(p - p'), (46) 

where 

^ = /pF^"'" s < 47) 

Inserting (f4*2"J) into eq. (jlPjl we find that 

Z(t) = J d n P A(p,p;t). (48) 

This form is analogous to the expression for Z(t) in terms of the Bloch density matrix 
(see, e. g. 0): 

C(x,x';t) = ^V:(x)^(x')e-* B " = (x^x'), (49) 

from which it is obvious that the quantity 

A(p,p';t) = {p\e~ t6 \p'), (50) 

whose trace is the partition function, is the counterpart of C(x, x'; t) in the momentum 
representation. 

As an example, we consider the case of n = 2, i. e. the Yang-Mills quantum 
mechanics (YMQM) with the potential 

V(x,y) = ^-x 2 y 2 , (51) 

which was the subject of numerous studies (see for a review). This potential leads 
to a non-separable and non-integrable Schrodinger equation. Its Fourier transform is 
given by 

a 2 h 4 d 2 d 2 

v «-rt = —MM S{p '- p) ' (52) 
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which lets (I46j) take the form 



J^4(p,p;t) + ^-A(p,p;t) 



g 2 h 4 d 2 d 2 



2 dp' 2 dp 



A(p,p , ;t) 



(53) 



p'=p 



From (p)3*j) one finds that the quantum corrections to the classical result for Z(t) generally 
begin at the order h for the YMQM model [Zj. In the same way, one may show that 
the quantum corrections to Z(t) for the harmonic oscillator begin at the order h 2 , in 
agreement with the result found by the expansion of the well-known exact expression 
for Z(t). The generalization of ()53j) to higher dimensions (n > 2) is straightforward. 



6. Conclusions 



In the present paper we have presented several different forms of the differential equation 
for the quantum corrections to the partition function Z(t). Some of these forms appear 
to be novel. Any one of these equivalent equations can be used to generate the gradient 
expansion of Z{t) for a wide class of potentials or could serve as the starting point for a 
numerical evaluation. The analogy between some forms of the equation and other well 
studied equations of theoretical physics, such as the Fokker-Planck equation, suggests 
the possibility of novel approaches or approximation schemes for the calculation of 
Z(t). As an example, we have used the integral representation of the Bloch equation 
in momentum space to give a simple proof of the fact that the quantum corrections to 
Z(t) for the anharmonic Yang-Mills quantum mechanics begins at the order h 4 We 
have also found a hitherto unknown correspondence with the Fokker-Planck equation 
occurring in supersymmetric quantum mechanics. 
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